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Abstract 

In this tutorial, I will first explain the need for finite calculus using an example sum 
I think is difficult to solve. Next, I will show where this sum actually occurs and why it 
is important. Following that, I will present all the mathematics behind finite calculus 
and a series of theorems to make it helpful before concluding with a set of examples to 
show that it really is useful. 
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1 How to Evaluate Y^x=i x 2 ? 

One of the problems we’ll 1 learn to address during this tutorial is how to mechanically 
(i.e. without much thinking) compute the closed form value of the sum 

X> 2 . (') 


While many students may already know the closed form answer to this particular 
question, our quest for the answer will lead us to techniques to easily evaluate nasty 
sums such as 


x=l y= 1 

(2) 

and 


X>2*- 

x=0 

(3) 

Since we’ll be using Yx=i x ' 2 as a motivating example, -v 
selves with the first few values of this function. 

ve’ll first familiarize our- 


n 

1 2 3 4 5 6 


1 5 14 30 55 91 


Now, we’ll try a few techniques to evaluate this sum before honing in on the correct 
answer. 

The astute student reading this tutorial may have surmised that we’ll make some 
connection with calculus at some point. Hence, let’s see what happens if we just 
“pretend” that this was an integral from calculus instead. Switching the Yh to a f 2 
sign gives 

^ x 2 = J x 2 dx. 

Using calculus, we can immediately solve the integral formulation. 



Sadly, trying n = 2 shows us that this is wrong. 

2 9 

E * 2 = |- 




Graphically, we can see what went wrong in this derivation. 

1 Since I like to keep things informal, if I say we in the text, I really mean you and I. 
interesting trivia: The integral symbol was actually a long letter S for “summa” [3] . So such a substi- 
tution makes sense at least alphabetically. 
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We want the area underneath the stepped line, not the area underneath the smooth 
curve. 

It is ironic that we can easily determine the area underneath the smooth curve 
using the infinite calculus but have trouble determining the much more simple area 
under the stepped line. In short, that is the goal of this tutorial - to derive a finite (or 
discrete) analogue of infinite calculus so the finite sum 

2 

x > 2 


is no more difficult to solve than the “infinite” sum 


x 2 dx. 


Besides finite calculus, another way to compute the value of x ' 2 i g to look it up 
in a book, or try and guess the answer. Searching through CRC Standard Mathematical 
Tables and Formuale [6], on p. 21 we find 


2 n(n + l)(2n + 1) 

= « ' 

I’m not sure I would have been able to guess that answer myself. 


2 The Computational Cost of Gaussian Elimi- 
nation 

Now we’ll venture off onto a bit of a tangent. Because I like to keep things applied, 
I’d like to show that our example sum occurs in practice. It is not an equation that 
popped out of the aether as an interesting mathematical identity, but rather emerges 
during a computational analysis of Gaussian Elimination - a fundamental algorithm 
in linear algebra. The Gaussian Elimination algorithm accepts an m x m matrix A 
and computes a lower triangular matrix L and an upper triangular matrix U such that 
A = LU. 
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Algorithm 1 Gaussian Elimination 
Require: A e M mXTn 
U — A, L — I 
for k = 1 to m — 1 do 
for j — k + 1 to m do 
Ijk = Mjk/Ukk 
Ujlrn = UjS'-m ~ ljkU k %m 
end for 
end for 


In this tutorial, I won’t give a full derivation of Gaussian Elimination, but simply 
give the algorithm (from [4]) and analyze its computational cost. 

The computational cost of an algorithm is the number of addition, subtraction, 
multiplication, and division operations performed during the algorithm. To analyze 
this cost, we’ll look at how much work is done at each step of Gaussian Elimination, 
that is, each iteration of the outer loop. In the case when k = 1, then we let j 
go between 2 and m. In the inner loop, we perform 1 division (lj k = Ujk/u k k), rn 
multiplications (IjkUk^k-.m), and nr subtractions. Since we run this inner loop rn — 1 
times, we have 

2 m(m — 1) + (rn — 1) - 2 m 2 - 1 < 2 m 2 

work on the first step. On the second step, we let j go between 3 and m. In the inner 
loop, we now perform 1 division, m — 1 multiplications, and m — 1 subtractions, for a 
total of 

2 (m — 2 )(m — 1) + (m — 2) < 2 (m — l) 2 . 

This pattern continues for each of the following outer steps. 

Thus, the total runtime of Gaussian Elimination is less than 

2 m 2 + 2(m - l) 2 + 2(m - 2) 2 + . . . + 1 = 2x 2 = 2 x 2 . 

X=1 X=1 

Since we looked up the answer to this sum in the previous section, we know that 
Gaussian Elimination takes less 

2g x 2 = 2 m ^ m + 1 > (2m + 

total additions, subtractions, multiplications, and divisions. Now, let’s start determin- 
ing how we can compute Y^=i x ' 2 ourselves. 

3 The Finite Calculus 

Before launching into the derivation of finite calculus, I should probably explain what 
it is. By analogy with infinite calculus, we seek a “closed form” expression for sums of 
the form 

£/(*)• 
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By closed form, we mean involving some set of fundamental operations, such as addi- 
tion, multiplication, exponentiation, and even factorials. Whereas in infinite calculus, 
we needed to compute the area under a function exactly, in finite calculus we want to 
compute the area under a sequence exactly. We call this finite calculus because each 
sum is composed of a finite (rather than infinite) set of terms. One way to think of 
finite calculus is calculus on the integers rather than the real numbers. 

In calculus, we used the notion of derivative and anti-derivative along with the 
fundamental theorem of calculus to write the closed form solution of 

J b f(x)dx = F{b)-F{a), 

where 

Our immediate goal for finite calculus (and this section) is to develop a fundamental 
theorem of finite calculus of a similar form. 


3.1 Discrete Derivative 


Before we can hope to find the fundamental theorem of finite calculus, our first goal is 
to develop a notion of “derivative.” Recall from calculus that 


= Ji 


f(x + h)~ f(x) 


Since we aren’t allowing real numbers, the closest we can get to 0 is 1 (that is, without 
actually getting to 0). Then, the discrete derivative is 


A/(z) = f(x+ 1) - f(x). 


(As an aid to those readers who are skimming this tutorial, we’ll repeat this and future 
definitions in a formal definition statement.) 

Definition (Discrete Derivative). The discrete derivative of f(x) is defined as 


A/(x) = f(x+ 1) - f(x). 


What can we do with our new derivative? Prom calculus, we found that it was 
simple to take the derivative of powers f(x) = x m . 


—x = mx 


Hopefully, we’ll find such a simple derivative for finite powers as well. 


Ax m = (x + l) m — x m ' mx m 1 . 


A quick check shows that 


Ax = (x + 1) — x = 1, 
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Unfortunately, this simple form does not work for x 2 . 

Ax 2 = (x + l) 2 — x 2 = 2x + 1 ^ 2x. 

Now, our derivate for x 2 was only off by 1. Is there any easy way we can fix this 
problem? If we rewrite our previous failed derivative as 

A (x 2 ) -l = 2x, 

we find ourselves must closer. Now, we can pull the —1 into the derivative by observing 
that A(x 2 — x) = 2x + 1 — 1 = 2x. We can factor x 2 — x as x(x — 1). Thus, we have 

A(x(x — 1)) = 2x. 

I’ll leave it to the reader to verify that 

A(x(x — l)(x — 2)) = 3x(x — 1). 

It appears we’ve found a new type of power for discrete derivatives! These powers 
are known as “falling powers.” 

Definition (Falling Power). The expression x to the m falling is denoted x — . 3 The 
value of 

x— = x(x — l)(x — 2) • • • (x — (m — 1)). 

Using the falling powers, we’ll now prove that falling powers are analogous to regular 
powers in finite calculus. 

Theorem 3.1. The discrete derivative of a falling power is the exponent times the 
next lowest falling power. That is, 


Ax— = 


Proof. The proof is just algebra. 

Ax— = {x + 1)— — x— 

= (x + l)a;(a: — 1) ■ ■ ■ (x — m + 2) — x(x — !)■■■ {x — m + 2)(x — m + 1). 

= (x + 1 — x + m — l)x(x — 1) • • • (x — m + 2) 

= mxttPA. 

□ 

Now, we’ll prove a few other useful theorems. 

Theorem 3.2. The discrete derivative of the sum of two functions is the sum of the 
discrete derivatives of the functions themselves. 

A (f(x) + g{x)) = A f(x) + A g{x) 

3 Note that other authors use the notation (x) rn to denote the falling power [5]. 


6 



Proof. The proof is straight from the definition. 

A (f(x) + g(x)) = f(x + 1) + g(x + 1) - f(x) - g(x). 

Rearranging terms, we get 

A (/(a;) + g(x)) = A/(x) + Ag(x). □ 

Theorem 3.3. The discrete derivative of a constant times a function is the constant 
times the discrete derivative of the function. 

A ( c /(®)) = cAf(x). 

Proof. Simply factor out the constant from the application of the definition of the 
discrete derivative. □ 


3.2 The Indefinite Sum and the Discrete Anti-Derivative 

With some handle on what our finite derivative does, let’s move on to discrete integra- 
tion. First, we’ll define some notation. 

Definition (Discrete Anti-Derivative). A function f(x) with the property that 
A f{x) = g{x) is called the discrete anti-derivative of g. We denote the class of func- 
tions satisfying this property as the indefinite sum of g(x). 

^2g(x)Sx = f{x) + C, 

where C is an arbitrary constant. g(x) Sx is also called the indefinite sum of g{ x). 

The discrete anti-derivative corresponds to the anti-derivative or indefinite integral 
from calculus. 

J g(x ) dx = f(x) + c ^ g(x) 5x = f(x) + C 

So far, all we’ve done is establish some notation. We have not shown any anti- 
derivatives yet. However, we can easily do just that. Recall that 

Ax m = mx^PA-. 


Using this fact, along with Theorems 3.2, 3.3, allows us to state that 


x— Sx = — — - + c. 

^ TO + 1 


Now, let’s work toward a fundamental theorem of finite calculus. First, we need to 
define the notion of a discrete definite integral or definite sum. 


Definition (Discrete Definite Integral). Let A f[x) = g(x). Then 


J2a(x)Sx = f{b ) - f{a). 
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With the discrete definite integral, the theorem we’d like to show (by analogy with 
infinite calculus) is 


^2g(x)dx = ^2g(x). 


Unfortunately, this isn’t true as a quick check shows. 


TxSx =- 1 =fiM = io 

i 2 |i 2 


But, )T)x=i x = 15, so the theorem isn’t true. However, we were very close. Note that 

5 4 

^ xSx = 10 = ^ x. 


This revised form gives us the correct fundamental theorem of finite calculus. 
Theorem 3.4. The fundamental theorem of finite calculus is 
b 6-1 

^2g(x)dx = J ^g(x )• 

Proof. Here we have more algebra for the proof. Let A f(x) = /( x + 1) — f(x) = g(x). 


J2 = J2 f ( x + !) _ /(*) 

= /(« + 1) - /(a) + /(a + 2) - /(a + 1) + . . . 

+ /(6)-/(6-l) 

= /W - /(a), 

after the telescoping sum collapses. 

Now f(b) — f(a) = Yla9( x ) by definition. □ 

3.3 Helpful Finite Calculus Theorems 

Now that we have our fundamental theorem, this section is just a collection of theorems 
to make finite calculus useful. The uninterested reader can skip to Table 1 at the end. 

One of the more useful functions from calculus is f(x) = e x . This special function 
has the property that 

D (e x ) = e x and J e x dx = e x + C. 

Our first theorem is that there is an analogous function in the finite calculus - a function 
whose derivative is itself. To find it, let’s “round” e. If we do this right, the analog of 
e should either be 2 or 3, right? Let’s see which one works. 

A (2 X ) = 2 X+1 - 2 X = 2 • 2 X - 2 X = (2 - 1)2* = 2 X . 

A (3*) = 3 X+1 - 3* = 3 • 3* - 3* = (3 - 1)3* = 2-3*. 

Thus, e corresponds to the integer 2. 



Theorem 3.5. The function 2 X satisfies 

A ( 2 X ) = 2 X and ^2 x 8x = 2 x + C. 

Let’s compute the general derivative of an exponent. 

A ( c x ) = c x+1 — c x = (c — 1 )c x . 

Because c is a constant in this expression, we can then immediately compute the anti- 
derivative as well. 

c x 8x = — b C. 

^ c- 1 

One of the important theorems in infinite calculus is the formula for -^(u(x)v(x)). 
Let’s find the corresponding formula for finite calculus. 

A (u(x)v(x)) = u(x + T)v{x + 1) — u(x)v(x) 

= u(x + l)v(x + 1) - u(x)v(x + 1) + u(x)v(x + 1) - u(x)v(x) 

= v(x + 1)A u(x) + u(x) Av(x). 

Now, we can use this derivative to write a discrete integration by parts formula. 

Theorem 3.6. 

u(x) Av(x) Sx = u{x)v{ x) — ^ v(x + l)Au(x) Sx. 

Proof. If we take the anti-derivative on both sides of 

A (u(x)v(x)) = u(x)Av(x) + v(x + l)Au(x) 

we get 

u{x)v(x) = ^ u{x)Av(x) Sx + ^ v{x + l)Au(x) Sx. 

Rearranging this equation yields the theorem. □ 

There are three facts I’d like to finish up with. First, let’s look at the derivative and 
anti-derivative of some combinatorial functions. Second, does our formula for falling 
powers and their derivative also work for negative powers? Finally, what about a;—? 

The binomial coefficients often arise in Combinatorics. Let’s compute A(^). If we 
use the well-known formula, (£) = Qlj) + (V) (fc) “ CV) = we can 

easily show that 

We can rewrite this as 



Since the falling power is only defined for positive exponent, you may have been 
wondering about negative falling powers. Let’s handle them now. We can derive a 
definition by looking at a pattern in the falling powers. 

x- = x(x — l)(x — 2), 
x- = x(x — 1), 

X- = X, 

afi=l 
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To move from x- to x- we divide by (x — 2). Then we divide by (x — 1) to go from x- 
to x-. When we continue this pattern, we get 


x— = x—/ (x + 2) = 


1 


(x + l)(x + 2) 


Definition (Negative Falling Power). 




{x 


1 

l)(x + 2) • • • (x + to) 


Let’s check our derivative rule for negative falling powers. 


(x + 2)(x + 3) • • • (x + m + 1) (x + l)(x + 2) • • • (x + m) 
(x + 1) - (x + m + 1) 

(x + l)(x + 2) • • • (x + m + 1) 

= rna: ~ m ~ 1 . 


Amazingly, I haven’t led you astray here. It does! This fact implies that our formula 
for the discrete anti-derivative for falling powers holds for negative falling powers as 
well. 

Well, our theorem holds, unless m = — 1. In that case, we’d get x ~ = if’ which 
is a problem. Recall from infinite calculus, we had to appeal to the log(x') function to 
integrate K There is no easy way of getting the correct function intuitively. 


Theorem 3.7. Let 


H x 



1 

x ’ 


for integer x. The function H x is the antiderivative of K 


Proof. 


A (H x ) = \ + 


1 

2 


1 11 1 _ 1 

x + 1 1 2 x x+l 


□ 


4 Making Finite Calculus Useful: 

Stirling and His Numbers 

If you’ve gotten this far, you’re probably thinking. “Great David, we’ve shown all these 
really neat facts about finite calculus. But, you still haven’t shown me how to solve 

i > 2 


as you promised way back in the introduction.” Well, let’s do that now! 
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f(x) A f(x) 


x— 

X— 


2 X 

c* 

a 

u(x) + v(x) 
u(x)v(x) 
u(x) Av(x) 


mi— 

0k 

2 X 

(c- !)<? 

LE) 

A u(x) + Av(x) 

u(x) Av(x) + v(x + 1)Am(x) 


E f( x ) Sx 

hT 

2 X 

(m+l) _ 

E u{x) Sx + E v(x ) Sx 
u(x)v(x) — E + 1) A u(x) Sx 


Table 1: Here we give a set a useful finite calculus theorems. These 
should make solving even nasty sums easy (or at least tractable). All 
of the discrete anti-derivatives have the constant omitted for brevity. 


In order to solve regular powers, we need to find some way to convert between x m 
and x— so we can use our integration theorems. 

Let’s see all the conversions we can come up with ourselves. 




x 1 = x- 
x 2 = X- + X- 
x 3 =??? 


So we find that the first few powers were easy to convert, but then the higher powers 
are not so obvious. Let’s see if we can work with 

x 3 = ax- + bx- + cx- 


to come up with a formula for a, b , and c. 

ax-+ bx- + cx- = [a(x)(x — l)(x — 2)] + [6(x)(a: — 1)] + [c(x)] 
= [ax 3 — Sax 2 + 2 ax\ + [bx 2 — bx] + [cx] 

= ax 3 + (b — 3 a)x 2 + (2a — b + c)x. 


If we want 


ax- + bx- + cx- = x 3 , 


then we need the following coefficients 


a = 1 (b — 3 a) = 0 (2a — b + c) = 0, 

or a = 1, b = 3, and c = 1. Hence, we have 

x 3 = x- + Sx- + x-. 

That was a lot of work just for x 3 . I’m not going to attempt x 4 , but feel free if 
you are so inclined. Instead, I’ll devote the rest of the section to proving the following 
theorem. 
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Theorem 4.1. We can convert between powers and falling powers using the following 
formula 



where {™} is a Stirling number of the second kind. 


Okay, that theorem was a lot to digest at once. Before we even start to prove it, 
let’s pick it apart. What the theorem is saying is that we can convert any power into 
a sum of falling powers as long as we can compute these special Stirling numbers (I’ll 
get to these next). So, if we wanted x A in terms of falling powers, this theorem tells us 
the answer: 


4 = {or + ur i+ \tr + {tr + \4 x 


Let’s address the Stirling numbers now. 


4.1 Stirling Numbers (of the Second Kind) 

If you haven’t seen them before, you are now probably wondering about these Stirling 
numbers. To start, I’ll just define them. 

Definition (Stirling Numbers). The Stirling numbers of the second kind, 4 denoted 

{:}• 

count the number of ways of partitioning n distinct objects into k non-empty sets. 

I’m guessing that didn’t help much and you are still scratching your head. Hopefully 
working with some of these numbers will help. First, let’s try and compute { ] } . This 
is the number of ways of partitioning 1 item into 1 non-empty set. I think it’s pretty 
clear that there is only one way to do this. So {*} = 1. What about {g}, or the number 
of ways to partition 1 item into 0 non-empty sets? Since each of the sets has to be 
non-empty, there are 0 ways of doing this, {J} = 0. We can generalize this argument 
and show that 



if n > 0. And what about {[}}? We want to divide 0 things into 0 non-empty sets. 
Shouldn’t this be 0 too? Strangly, no, {jj} = 1. If you have trouble accepting this fact 
logically, just consider it a definition. 

Another set of simple Stirling numbers are {”}. These are the number of ways 
of dividing n objects into n non-empty sets. Since the number of objects equals the 
number of sets, there can only be one object in each set, so {”} = 1. 

Now, let’s look at a more complicated Stirling number, { 2 }. Suppose our three 
items are 0, and < v ) . When we divide these three card suits into 2 non-empty sets 

4 There is another set of numbers called Stirling numbers of the first kind which are slightly different. 

See [5] for more information. 
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we get the following sets. 


{*,0},{9} 

{*,9},{0} 

{0,9}, {4k} 

So { 2 } = 3, since there are three ways to put 3 items into 2 non-empty sets. 

We still don’t really know much about Stirling numbers though. Let’s prove a 
theorem to help us general new Stirling numbers from previous ones. 



Proof. We’ll give a combinatorial proof of this identity. In a combinatorial proof, we 
choose a particular quantity and give two ways of counting this quantity. Since these 
two different ways count the same thing, they must be equal. In this proof, we’ll count 
the number of ways to partition n distinct objects into k non-empty sets. 

The obvious way to count this quantity is by using our previous definition of Stirling 
numbers. By this definition, {£} counts the number of ways to partition n distinct 
objects into k non-empty sets. 

Now, let’s count this quantity using a different way. I like to call this the “weirdo” 
method. Let’s pick out one of our distinct objects and call it x, the “weirdo.” We can 
now count the number of ways to partition n distinct objects into k non-empty sets by 
keeping track of where the “weirdo” goes. The first possibility is that the “weirdo” is 
in a set all by itself. In this case, there are {^Ii} ways for the other n — 1 objects to 
be partitioned into k — 1 non-empy sets. 

If the “weirdo” is not in a set by itself, then it must be inside another non-empty 
set. So the others must have been divided into k non-empty sets. There are {"jj 1 } 
ways of partitioning them. For each of these partitions, we can insert the “weirdo” into 
any of the k sets and get another partition. So there are k{ n f, 1 } total ways putting the 
“weirdo” into these existing sets. Since we have completely dealt with the “weirdo”, 
there are {^Ii} + k{ n ~f ways to partition n distinct objects into k non-empty sets. 

Thus, 

a-LJM"'} 

This proof is a little challenging. Let me just go through it with our 3 items from 
before. Let 4k be the “weirdo.” (It does look a little strange, doesn’t it?). If 4k is in a 
set by itself, then we need to partition •Os’v’ between the 1 remaining set. There is no 
way to divide them between one set, so there is just one way of putting them in a set 
together. We get the following partition 

{4k}. 

The other case was that 4k was not by itself. So we partition <C>, U 1 into 2 sets. 
There’s still just one way of doing this {0}, {9}. We can insert 4k into either of the 
sets in this partition, so we get 2 more ways of partitioning the items into two sets. 
{*,0},{9} 

{*,9},{0} 

Using this theorem, we can build the Stirling number triangle (like Pascal’s trian- 
gle). See Figure 1. 
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1 


0 1 
0 1 1 
0 13 1 

0 17 6 1 

0 1 15 25 10 1 

0 1 31 90 65 15 1 

Figure 1: A table of Stirling Numbers of the second kind {”} which 
induces the Stirling triangle. The value of {”} is the k + 1th number on 
the n + 1th row, e.g. { 2 } = 7. To build the triangle, use Theorem 4.2. 
This theorem says that the value in a particular position is k times the 
number up and to the right, plus the number up and to the left, e.g. 
7 = 2*3 + 1. 


4.2 Proving the Theorem 

In order to prove Theorem 4.1, we’ll first prove a helpful lemma, then do a lot of 
algebra. 

Lemma. A helpful lemma 

x . x k = x k+A + kx k 


Proof. 


= (x — k) ■ x- + k ■ x- 
- x ktl + k-x k . 


□ 


Now, recall Theorem 4.1. 



Proof. We’ll prove this by induction on n. 

The base case, n = 1 is trivial, x 1 = + {l}®- = x. 

Our induction hypothesis is that x n = Y2k = 0 holds for all n < r. 
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(See the end for a discussion of the steps in the following work.) 



In the first step, we applied our induction hypothesis, then we moved the x inside 
the sum. Following that, we applied the lemma we just proved. We then split the sum 
into one sum based on k + 1 and the other based on k. Then, we pulled out the {£} = 1 
and {q} = 0 terms from the sum. When we put things back together in the sum, we 
use Theorem 4.2. □ 

Excellent! Now we have a means to convert between x m and falling powers! 

4.3 Computing Stirling Numbers 

We can derive a formula to compute Stirling numbers without the recurrence relation 
derived in Theorem 4.2. This derivation is somewhat complex, and I’d encourage 
readers to skip it unless they are really interested. The important part of this section 
is the statement of the theorem. 

Theorem 4.3. 



Proof. We’ll prove this using a combinatorial argument. The objects we are count- 
ing are the number of surjections between an n-set X and a k - set Y . Let Y = 
{j/i, ?/ 2 , • • • , Uk}- Recall that a surjective mapping is a mapping onto the set Y, so 
some element in X must map to each element in Y . 

First, we can count this using /c! {))} . To see this fact, consider that after the 
partitioning of n into k non-empty subsets, we have a partition such as the following 
one 

{xi,x 5 ,x 9 }, {x 2 ,X±}, . . . 


1 ! 



If we apply an ordering to the sets comprising these partitions, i.e. 7 Ti, 7T2, . . . ,7Tfc. then 
we can construct a surjective map from X to Y using the map 

/(*) = Vj | X G 7 Tj, 

that is, the set 7 Tj is the inverse image of yj. Since there are k\ orderings of the numbers 
between 1 and k, we get that there are surjections. 

We can also count this using an inclusion-exclusion argument. If we take the total 
number of mappings from X to y, we want to subtract all the mappings that do not 
map to at least one element in y. We can do this by look at all maps that exclude 
a particular element y*. To do this, we take the total number of maps between X 
and Y and subtract the number of maps that exclude an element y t . However, now 
we’ve subtracted some maps twice (maps that exclude y, and y 3 ) , and we then need 
to add them back. So we then add the number of maps that exclude two elements, 
then subtract the maps that exclude at least three (for a similar reason), etc. This 
strange sum actually counts all of our elements. For more information about a powerful 
generalization of this proof technique, see [5] . 

There are (^) ways to pick the i elements from k to exclude and there are ( k — i) n 
maps that exclude one set of i elements from Y . Thus, we get the number of surjective 
maps from X to Y is 

(;)(*-*)" 

We divide both sides by A:! to complete the proof. □ 


all maps exclude 1 exclude 2 


5 Examples Galore... 

In this section, we just provide a few in-depth examples to show the utility of finite 
calculus and how it really takes the thought out of complicated sums. 

5.1 EEi* 2 

Since I still haven’t gotten around to showing how to derive the closed form solution 
to this sum, let’s do that now. It’s a nice quick warm-up to the following examples. 

n n + 1 n+1 

x 2 = ^ x 2 8x = ^ x- + x-Sx = x-/3 + x-/2|” +1 = (n + l)-/3 + (n + 1)— /2. 

5.2 A Double Sum 

As promised in the introduction, let’s compute Ylx=i J2y=i( x +y) 2 using finite calculus. 
First, we’ll convert to falling powers and take the anti-derivative, then evaluate the 
inner sum and repeat the process. 
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n n n+1 n+ 1 n+1 n+1 

= J2Jl( x+ y') 26 y 6x = J2Jl( x + y') 1+ ( x+ y') l6 y 6x 

X=1 y = 1 11 11 

n+1 

= ^(cc + n + l)-/3 + (x + n + l)-/2 — (x + l)-/3 — (x + l)-/2 Ac 

= (2 n + 2)-/12 + (2 n + 2)^/6 - (re 4 - 2)^/6 - (n + 2)^/3 

Oops. While that was correct, we skipped one step. What is the J2( x + c)— 8x1 
We assumed it was the same as Yl x ~^ x without proof. Looking at the proof of 
Theorem 3.1, it is true since we just add a constant to the falling power. 

5.3 Average Codeword Length 

(This example, without the story, comes courtesy of CS 141 at Harvey Mudd Col- 
lege [2]). Suppose you are the director of programming for one of the Mars rovers. 
These rovers accept command codewords of up to 10 bits. A brief summary of the 
commands are listed below. 

0 stop 

1 take photo 

00 move forward 

01 move backward 

10 move left 

11 move right 

001101 take picture of green alien 
1111111111 self-destruct 

There are a few interesting aspects to this table. First, 0 and 00 are different com- 
mands! Second, all possible commands are used and there is no redundancy between 
the command codes. 

After you have designed and implemented this codeword system, your boss calls 
you on the phone and wants to know the average codeword length. Since NASA pays 
interplanetary communications costs per bit, he really wants to know how much money 
your codeword system will cost. 

While still talking to your boss, you quickly jot down a few notes. 

First, for a codeword of exactly x bits, you note that there are 2 X possible codes. 
Hence, the total number of codes is 



Second, the total length of all codewords of exactly x bits is x2 x . Aha! Now you 
can easily compute the average codeword length. 

Ex=i x2X 
Ex-i 2 * ' 
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While you’ve already solved the denominator, the numerator is a little more chal- 
lenging. However, since you have (obviously) read all portions of this tutorial you know 
about finite calculus and, unabashed, proceed ahead. 

10 11 

y] x2 x = y x2 x Sx. 

Now, you remember this little theorem about “discrete integration by parts” (The- 
orem 3.6). 


y x2 x Sx = x2 x - y 2 X+1 Sx = x2 x - 2 X+1 = 2 x (x - 2). 

Then, 

10 

y®2*- 2 x (x — 2)| J 1 = 9 • 2 11 + 2. 

The average codeword length is then 
q . 9 11 _i_ o 

2ll _ 2 = 9.00978 « ». 

After a brief pause in the conversation with your boss, you definitively state: “Sir, 
the average codeword length is 9 bits to two decimal places.” 


6 Conclusion 

I hope you’ve enjoyed learning about finite calculus as much as I’ve enjoyed writing 
about it. I think that finite calculus is a very nice set of mathematics that is quite 
easy to understand and apply. Hopefully, you now understand it and can apply finite 
calculus to your own work if the need arises. If you are still interested in the subject, 
please see the excellent book, Concrete Mathematics [1], by Ronald Graham, Don 
Knuth, and Oren Patashnik. In that book, the authors provide a similar derivation of 
finite calculus and show many other cases where it helps solve yet another nasty sum. 
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